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Abstract 

Hypothesis elimination is a special case of Bayesian updating, where each piece 
of new data rules out a set of prior hypotheses. We describe how to use Grover's 
algorithm to perform hypothesis elimination for a class of probability distribu- 
tions encoded on a register of qubits, and establish a lower bound on the required 
computational resources. 



1 Introduction 

In the standard setting of Bayesian inference one starts from a set of hypotheses HI = {h] 
and a set of possible data D = {d]. Hypotheses and data are connected via conditional 
probabilities p{d\h), known as the model. Given a prior probability distribution p{h), the 
updated probability p{h\d) of the hypothesis h to be true given that the piece of data d 
was observed is obtained from Bayes's rule, [IJ 

P(m = ^"'"''Tlt'^M ■ (1) 

In this paper, we consider the problem of hypothesis elimination, which is a special case 
of Bayesian updating where the model is of the form 

(ri\h\ — i ^ if /i is ruled out by d , , . 

\ Ch otherwise. 

The positive constant Ch does not depend on d and is determined by normalization. We 
assume that there is a finite number, N, of hypotheses, which we label to — 1, i.e.. 
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EI = {0, 1, . . . , — 1}. Furthermore, we assume that the prior p{h) has been obtained 
by hypothesis ehmination from an initial uniform prior distribution on H. 

There are several possible ways of encoding a probability distribution p{h) on a 
quantum register. Here we represent p{h) by the state 



N-l 



h=0 

where | h) are computational basis states j2] of a register formed of [log2 N] qubits. 

To be specific, we assume that the prior p{h) has been obtained by A; — 1 hypothesis 
elimination steps {k > 1). We thus assume that the prior is given in the form of a 
sequence of oracles oi, . . . , Ok-i, where 

X _ J if /i is ruled by the data in step j , , s 

\ 1 otherwise. 

Likewise, we assume that the model Eq. Q is given as an oracle Ok, where Ok{h) = 
c'^^p{d\h). The sequence of oracles Oj gives rise to a new sequence Oi, ... ,0k defined by 

For each Oj, we define the set of solutions, Qj = {h : Oj{h) = 1}, and denote by Mj the 
corresponding number of solutions, 

M, = Y,0,{h). (6) 

For j G {/c — 1, fc}, we now define the states 



The state I'^k-i) encodes the prior p{h), and encodes the posterior p{h\d). 



2 Quantum hypothesis elimination 

The problem of hypothesis elimination now takes the following form: Given a quantum 
register in the prior state {"^k-i) and given the oracles Oi, . . . , Ok, transform the register 
state into the posterior state I'^k)- One can use Grover's algorithm j3j to solve this 
problem as follows. 

Define a quantum oracle corresponding to Ok via 

dk\h) = i-if^^'^lh) . (8) 

There are standard techniques [2^ to implement Ok in the form of a quantum circuit. 
The Grover operator G{Ok) associated with the oracle is then defined as 

G'(Ofe) = (2|^o)(^o|-l)0fc, (9) 
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where 1 is the identity operator and 

1 

w^^Ew (10) 

is the equal superposition state. The posterior state can now be prepared by 
repeated apphcation of the Grover operator G{Ok) to the equal superposition state I^E'o)- 
This requires [tt / A)^J N / calls of the oracle Ok- 

Notice that the hypothesis elimination algorithm outlined above makes no direct use 
of the prior state \^k-i)- This raises the following question: Is it possible to reduce 
the number of Grover iterations (and therefore oracle calls) required to prepare \^k) 
by starting from the prior state |^'fc-i) instead of the equal superposition state l^'o)? 
In other words, can one make use of the computational effort that went into preparing 
the prior state I^E'fc-i) in order to obtain the posterior state \^k) more efficiently? As 
already suggested by the results in Ref. j4j, the answer to this question is negative. Here 
we prove the following result. 

Consider the family of oracles O that consists of Ok-, Ok-i and all possible combina- 
tions of Ok and Ok-i (see Eq fj29|) for the precise definition of this family). Now consider 
all possible algorithms that consist of applying the corresponding Grover operators 

Aabc... = --- G{0,)G{0,)G{Oa) , Oa, Ob, 0„ ■ ■ ■ G O . (11) 

Then, in the limit of large A^, any algorithm Aabc... requires at least {V2/8)^N/Mk 
oracle calls from the above family to convert |^E'fc_i) into {"^k)- In other words, making 
direct use of the prior state I^E'fc-i) does not improve the asymptotic cost of 0{^y N/Mk) 
oracles calls required to prepare l^^k)- 



3 Proof 

Consider an oracle O which accepts M out of the total hypotheses h: 

N-l 

Y,0{h)=M. (12) 

We shall call such hypotheses good, as opposed to bad hypotheses that are rejected by 
the oracle. Using different notation ^ for the amplitudes of good and bad hypotheses, 
we have that after t consecutive applications of the Grover operator G{0) an arbitrary 
quantum state 

1*^'^')= E 9T\h)+Y.^^'\h) (13) 

good h bad h 

is transformed into 

good h bad h 



3 



Let (7'" and 6™' be the averages of the initial amphtudes corresponding to the good and 
the bad hypotheses respectively: 

9'"' = iEfi"'. s'-'-F^E'r. (15) 



h bad h 

and similarly for the final amplitudes 

f-^sEs?. 5«" = ^^6r, (16) 

good h bad h 

Let US also define 

^gr = gr - , = k - ■ (17) 

In other words, A(yf™' and define the features of the initial amplitude functions g^™ 
and relative to their averages (7™' and Biham et. al. |1] have shown that the 
change of the amplitudes is essentially determined by the change of the averages: 

^fin _ h^^' + {-ly M)f , (18) 

where the averages g^^ and W'^ are given as follows. Define 

f 2M\ , , 

uj = arccos II — j , (19) 



a = ^/|5ini|2+|^ini|2_^^ (20) 



arctan^W— — . (21) 

I 51111 y N - M ' ^ ^ 



The averages are given by 



N-M 

^ = y " sm(^^ + 0), 

b^^ = a cos{iut + (f)) . (22) 
Let us also define the separation of the averages 

^ini ^ -ini _ pni ^ ^fin ^ -fin _ gfin _ (23) 

Equations (fTH|) imply that after t applications of the Grover operator the values of indi- 
vidual amplitudes gf^ of accepted hypotheses do not change relative to their average g^'^ 

- = AC = gf - . (24) 

The same is true for the rejected hypotheses if A6"' = or t is even. This observation 
suggests that the action of any algorithm of the type 

Aaabbcc... = ... G\0,)G\0,)G\Oa) (25) 



4 



can be analyzed by looking at the changes of the average amplitudes of the accepted and 
rejected hypotheses relative to various oracles Oq, Ob, .... Before we proceed with this 
analysis let us first clarify what kind of oracles are relevant to this problem. 
Let / be a real-valued step function which takes only three values, 

if Ok{h) = l, 

f{h)={ h if Ofc_i(/i) = l, butOfe(/i)=0, (26) 
if Ofc_i(/i) = 0. 

Let us also require that Ylih=o — 1; introduce a quantum state |^'(/)) that is 

defined by / in a natural way, 

l^(/)) = E/(^)l^)- (27) 

h=0 

Evidently, both |\E'fc) and {"^k-i) can be written in this way. It follows from the above 
discussion that the action of the operator G^{Ok) on the state (P7|) can be completely 
described by the changes of /i, /2 and f^. Moreover G'^iOk) preserves the value of 
^k-i = f2 — /a- Similarly, the action of G'^{Ok-i) on (jTTjl can be completely described 
by the changes of /i, /2 and /a, and it preserves the value of Sk = fi — f2- In general, 
for any oracle O, the corresponding operator G^(0) preserves the amplitude differences 
between any two hypotheses for as long as either both hypotheses are accepted or both 
are rejected by O. However, if /good and /bad denote the amplitudes of an accepted and 
a rejected hypothesis, respectively, then the difference 6 = /good — /bad is changed by an 
amount A which satisfies the inequality 

|A| < . (28) 

where M is the number of accepted hypotheses with respect to O (see the Appendix). 
Using this inequality the action of any algorithm of the type ()25|) can be analyzed by 
calculating how individual changes of 8^ and accumulate during the action of the 
algorithm. In order to convert I^E^fc-i) into \^k) the net result of such changes must be 
sufficient to increase 5^ from to l/i/Mfc and decrease from \j ^fM]~[ to 0. 

It follows that all oracles that are relevant for this task can be obtained from and 
Since each oracle is completely characterized by the set of acceptable hypotheses, 
the relevant family of oracles generated by and Oh-\ can be written out as oracles 
that correspond to the sets 

ilfc, ^fc-i7 ^fe-i n , and the complementary fifc, fi^.i , fi^ U ■ (29) 



Let us consider the first three oracles from the family defined by the sets ()29|). namely 
the oracles that accept hypotheses from the sets 

^fc, ^fc-i7 ^fc-i n ilfc . (30) 
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Oracles that correspond to the complementary sets in ()29|) can be analyzed in a com- 
pletely analogous manner. The oracles corresponding to the sets flHUI) are Ok, Ok~-i and 
0^_^, where 

OlAh) = {\ (31) 
U otherwise . 



Using the inequality (|2H|) we see that, regardless of its position in the algorithm, 
the operator &{0k) changes 5k by at most di = 4:^/2 /y/N — Mk without changing Sk-i- 
Similarly, the operator ^(Ok-i) changes Sk-i by at most c/2 = 4v^/ ^^A^ — Mk-i without 
affecting the value of 6k. The Grover operator with the combined oracle G{0^_^) changes 



both 6k and 6k-i by an equal amount that does not exceed = 4\/2/ y — {Mk-i — Mk). 

Let N{0) be the number of times that the oracle O is called within the algorithm. 
We would like to find a lower bound on the total number of oracle calls, A'totai = N{Ok) + 
N{Pk~i) + N{p^_^), that is needed by the algorithm to convert l^&fc^i) into \^k). Let 
X, y and z be the number of times that the operators G'^(Ofc), G'^{Ok-i) and G'^{0'l_^ 
appear in the algorithm. Then A^totai is bounded from below by the minimal value of 
2{x + y + z) subject to the constraints 

xdi + zd'i > U , 

yd2 + zd3 > V, (32) 

where U = l/i/Mfc and V = Xj \JMk^\ are the required changes of 6k and 6k-\ respec- 
tively. This is a simple linear optimization problem that should be considered in the 
nonnegative x,y,z octant. Keeping Mk and Mk-i constant we obtain that in the limit 
of large the value of A^totai is approaching [yf^j N j Mk . 

For any algorithm Aahc... the action of the corresponding sequence of Grover operators 
can be rewritten in the form similar to that in Eq. Using the definition of the Grover 
operator, we have 

G(o.j . . . G(o,3)G(o,jG(o,j = (i6„j . . . (ia3)(iaj(iaj , (33) 

where A = 2|\E'o)(\I^o| - 1- Using the fact that {Okf = 1 and = 1 we obtain 

{AdaJ...{Ada,){Ada,){Adj = Ada^...Ada,AOaMOaJ'\Ada,Y 

= Ad,„...Ad,,A{da,a,A){Ada,)\ (34) 

where Oajaa = OajOa^. Denoting O 

= Ili=i we proceed 

i6,„ . . . Ada,AiOa,a,A)iAda,f = ... Ada,AiOa,a,A)-\da,a,A)\Ada,y 

){Oa,a,Af{AO 
= RG (Oa^.„a„) . . . G (Oaiazaa)^ (Oaiaa)^' (OaJ , 

(35) 

where the -|- and — signs are chosen for odd and even values of n, respectively, Oai...aj 

denote classical oracles that correspond to the quantum oracle Oai...aj , and /2 is a residual 
operator 

B = I ^('^'^i-'^n) if n is even, , . 

I G-\Oa,...a,J Otherwise. ^ ^ 
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Since all oracles, Oaj...a., belong to the family associated with the sets (1221), the above 
arguments allow us to derive a bound on the minimum number of oracle calls from 
this family that are required to convert {"^k-i) to I^I/fc). Indeed, transformations between 
Eqs. (jH^ and (jHUj) at most double the number of oracle calls that are used by the original 
algorithm. To be more precise, if n is the number of oracle calls used by the original 
algorithm (see the left-hand side of Eq. (jHSl))) then the equivalent modified algorithm, 
defined by the right-hand side of Eq. (j35j) . requires at most 2n+ 1 oracle calls. It remains 
to note that after the application of the residual operator that concludes the modified 
algorithm one has to arrive at the target state or, which is equivalent, the algorithm 

G'±'(0.,...J . ..G'{0,,a,a,)G-'{0,,a,)G\Oa,) (37) 

must prepare the state \ipk) = G^^{Oai...a„)\'^k)- In the limit of large the state 
l^pk) coincides with l^^k)- Using an analysis analogous to that of algorithm ^Hjl we 
therefore conclude that, in the limit of large A^, the algorithm (jHTf) requires at least 
(\/2/4)^iV/Mfc oracle calls to convert '^k-i into iV'fc)- The original algorithm, therefore, 
will need, asymptotically, at least {y/2/8)^/N/ Mj. oracle calls to convert \'^k-i) into 



Appendix 

Using the notation of Eqs. 



, the inequality (j28|) can be written as 



VN-M 

where the number of iterations is t = 2. By definition, we have 



l^fin _ ^inil 



a 




M_ 

N 



sin(2co' + ^) — sin,^ 




ttA/— I 2 sintu cos(co' + ,^)| 



V M 



sm — 
2 



Since < uj < n 



and therefore 



CO 



sm 



1 — COSCiJ 




, fin A ini 



An < 4a. 



(38) 



(39) 



(40) 



(41) 



A bound on a can be obtained using the fact that, for any ai, 02, . . . , ajv/ G M such that 
= 1) we have 

N n.i.] = — 

M 



A;=l 



M 



(42) 
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This can be easily shown using the method of Lagrange multiphers. Using ()42j) we obtain 
max q = , , and max b = = . (43) 

It then follows that 



/I U1 ,1 -12 ^ (AA\ 

a < W maxor + maxfifr— — = —= (44) 

Combining this bound with Eq. pT|l we obtain Eq. ()38|1 as intended. 
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